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ABSTRACT 

NASA Langley Research Center is investigating a guided wave-based defect detection technique 
for as-fabricated carbon fiber reinforced polymer (CFRP) composites. This technique will be 
extended to perform in-process cure monitoring, defect detection and size determination, and 
ultimately a closed-loop process control to maximize composite part quality and consistency. The 
overall objective of this work is to determine the capability and limitations of the proposed defect 
detection technique, as well as the number and types of sensors needed to identify the size, type, 
and location of the predominant types of manufacturing defects associated with laminate layup 
and cure. This includes, porosity, gaps, overlaps, through-the-thickness fiber waviness, and in- 
plane fiber waviness. The present study focuses on detection of the porosity formed from 
variations in the matrix curing process, and on local overlaps intentionally introduced during layup 
of the prepreg. 

By terminating the cure cycle prematurely, three 24-ply unidirectional composite panels were 
manufactured such that each subsequent panel had a higher final degree of cure, and lower level 
of porosity. It was demonstrated that the group velocity, normal to the fiber direction, of a guided 
wave mode increased by 5.52 percent from the first panel to the second panel and 1.26 percent 
from the second panel to the third panel. Therefore, group velocity was utilized as a metric for 
degree of cure and porosity measurements. 

A fully non-contact guided wave hybrid system composed of an air-coupled transducer and a laser 
Doppler vibrometer (LDV) was used for the detection and size determination of an overlap. By 
transforming the plate response from the time-space domain to the frequency-wavenumber 
domain, the total wavefield was then separated into the incident and backscatter waves. The 
overlap region was accurately imaged by using a zero-lag cross-correlation (ZLCC) imaging 
condition, implying the incident and backscattered waves are in phase over the overlap boundaries. 

1. INTRODUCTION 

The NASA Advanced Composites Project seeks to develop and transition technology that will 
reduce the timeline required for development and certification of new aircraft structure that utilizes 
advanced composite materials. One area of concern is the formation of porosity and fiber waviness 
defects in composite structural parts during various processing stages, from incoming material 
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handling to the completion of laminate curing in an autoclave/vacuum press. Defects have 
detrimental effects on the quality of the structural parts, often leading to higher cost and significant 
delay in parts production and utilization. 

In the composite industry, cure cycles are developed from a “trial and error” or a more effective 
“processing science” approach to reduce the final porosity level in the composite laminate [1]. 
Several studies have been published on the effects of porosity and their detection using non- 
destructive inspection techniques. Porosity is defined as “large number of microvoids each of 
which is too small to be of structural significance or to be detected individually by a realistic 
inspection technique, but which collectively may reduce the mechanical properties of the 
components to an unacceptable degree” [2]. Typically, porosity occurs because of entrapped air, 
moisture, or volatile products during the curing cycle [2]. In a special case with pure resins only, 
ultrasonic velocity has been used to infer the degree of cure because of its effect on the modulus 
of the resin [3-5]. In CFRP, it has been shown that as the porosity level increases, the ultrasonic 
velocity of the wave propagating perpendicular to the carbon fibers decreases [6]. Ultrasound can 
monitor end of cure based on when the time delay between the pulse and echo of the transducer 
plateaus (i.e. the ultrasonic velocity ceases to increase) in graphite/epoxy composites, [7] and 
epoxy matrices [8]. Non-contact, air-coupled transducers can be used for ultrasonic velocity 
measurement in resins when line of sight is available during curing [9]. Using mounted 
piezoceramic actuators and sensors can alleviate the unreliable coupling between an ultrasound 
transducer and the mold which can occur, especially during the phases of heating and cooling [10]. 
Other ultrasonic phenomena have also been used for monitoring degree of cure including 
attenuation [2,5-7, 11], instantaneous phase, and the “mean value of each frequency curve weighted 
by the maximum corresponding spectral amplitude” [12]. 

Based on the premise that the level of porosity may be dictated by the change of the ultrasonic 
wave speed, this work examined the group velocity of guided waves in several panels with 
different final degrees of cure and lower level of porosity. The significance of this approach was 
that unlike conventional ultrasound that provides information about the quality of the part directly 
underneath or nearby the ultrasonic transducer in a point-by-point (discrete) inspection; the guided 
wave approach interrogated a continuous wave path through the thickness of the panel along the 
line from actuator to sensor. In addition, the experiment was designed such that it may be 
transitioned into a real-time detection system during cure. If the level of porosity can be estimated 
during curing, the cure cycle can be dynamically adjusted based on the measurements. This would 
prevent the operator from having to follow a prescribed, or non-optimized cure cycle. 

Ply gaps and overlaps occur primarily at tow boundaries during the layup process and can lead to 
surface waviness [13]. Their effect on the ultimate strengths of composite laminates have been 
experimentally investigated [14]. Several optical imaging techniques have been or are being 
developed for the detection of gaps and overlaps during layup from an automated fiber placement 
(AFP) machine. The detection system can be manual or fully automated [13,15,16]. However, 
each of these systems cannot detect gaps and overlaps in a panel that has already been 
manufactured. Eddy currents have been used in a controlled research environment for the 
detection of gaps and overlaps on uncured prepreg [17]. 

For this study, the local defect detection and size determination for the overlap was devised by a 
fully non-contact laser ultrasound system based on a guided wave-based approach. The system 
was composed of a stationary air-coupled transducer, operating in the vicinity of the overlap 
region, and a laser Doppler vibrometer (LDV) that captured the wavefield by scanning the region. 



Here, the total signal energy was first mapped in the region to show the area that was primarily 
illuminated by the air-coupled transducer. The wavefield, in time segments, in the illuminated 
region was then transformed into the frequency- wavenumber domain. In this domain, the incident 
and backscatter wavefield can be readily separated. Frequency-wavenumber filtering has 
previously been used for improved damage imaging of longitudinal slits in aluminum plates [18- 
20], disbonding in a composite wing [21], and delamination in a multi-layer composite plate [21]. 
The imaging metrics utilized were root mean square distribution of the reflected signal [18], 
cumulative standing wave energy [19], zero-lag cross-correlation (ZLCC) [20], standing wave 
filter [21], and Laplacian image filter [21]. For this work, the ZLCC imaging condition was 
utilized to reveal and image the overlap region. 

2. EXPERIMENTATION 


2.1 Porosity Detection 

Three panels of twenty-four ply unidirectional Hexcel® IM7/8552 prepreg were fabricated in a 
vacuum press according to the cure cycles shown in Figure 1. The panels were 305 mm x 305 mm 
with nominal thickness of 3 mm. Each subsequent panel had a higher degree of cure than the 
previous one, based on the temperature and pressure experienced during the cure, also shown in 
Figure 1. 



Cure Time (mins) 

Figure 1. Cure cycle for Panel #1, #2, and #3. 

After cure, each of the three panels were C-scanned (Figure 10) with a Sonix, Inc. immersion 
ultrasound system in pulse-echo mode using a Panametrics 1 5 MHz transducer. After C-scan, the 
panels were investigated by the guided wave approach. A guided wave was excited into the 
laminate using a five-peaked toneburst signal with a center frequency of 300 kHz emitted from a 
waveform generator (Agilent Technologies: 81150A) to an amplifier (Krohn-Hite Corporation: 
Model 7602M) to a piezoelectric disk (Steiner and Martins Inc.: SMD07T02R4I2WL). The 
laminate response was recorded by another piezoelectric disk in a pitch-catch configuration 
(Figure 2). To measure a wave propagating normal to the fiber direction (i.e., 90 degree direction), 
shown in Figure 3, each of the piezoelectric disks in the bottom row served as an actuator. The 
laminate response was recorded by the three piezoelectric disks along the same vertical column. 
This process was repeated for each of the four colu mn s giving a total of twelve measurements. 
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Figure 2. Experimental setup for deteeting Figure 3. Aetuation and sensing locations 
porosity level related to the ultrasonic wave normal to the fiber (ninety degree) direction, 
velocity of a 24-ply unidirectional composite 
using piezolelectric disks. 


2.2 Overlap Detection 

A panel was laid up by hand using Hexcel® 1M7/8552 unidirectional prepreg. The layup was 
[O 4 / 9 O/ 9 O 2 / 9 O/O 4 ] where in the middle two layers overlap and gap defects were intentionally 
introduced. The entire panel was 150.4 mm x 150.4 mm x 1.524 mm (nominal) (6 in. x 6 in. x 
0.06 in. (nominal)) and the overlap and the gap were of dimensions 63.5 mm x 6.35 mm x 0.25 
mm (2.5 in x 0.25 in. x 0.01 in.) (Figure 4). The panel was cured in a vacuum press following the 
cure cycle recommended by Hexcel®. The temperature and pressure scheme for the cure cycle is 
similar to that of Panel #3 in Figure 1. 
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Figure 4. Layup of panel with intentionally introduced overlap and gap defects near the neutral 

ply surface of the composite laminate. 

An unfocused circular air-coupled transducer with an active diameter of 25 mm (Ultran Group: 
NCG200-D25) was used to excite a five-peaked toneburst acoustic wave with a center frequency 
of 250 kHz (Figure 5). The excitation area on the plate was 26.5 mm x 25 mm. The location of 
the air-coupled transducer was stationary while the out-of-plane velocity of the Ao Lamb wave was 



recorded in a 1 mm x i mm grid inside the scan area (46 mm x 73 mm) (Figure 6 ) by a LDV 
(Polytec Inc.: OSV 500) (Figure 5). 



Figure 5. Experimental setup for detecting 
overlap region. 


Figure 6 . Excitation area from air-coupled 
transducer and LDV scan area (in blue) and 
imaged area (outlined by red rectangle). 


3. RESULTS AND DISCUSSION 

3.1 Porosity Detection 

The group velocity of the guided wave, c^, was calculated as the distance between the actuator 
and the sensor, , divided by the time of flight (TOF) (Eq. (1)) 


where is the time of actuation and is the time of arrival at the sensor. 

To determine and , the following algorithm was implemented. First, the time at which there 
is minimum response was determined and recorded as • Then, the time at the maximum peak 

prior to was determined and recorded as I . The reason must be prior to 

A «t‘'l|min “VI Imax “VI Ij^^x ^ “VI|min 

is that the asymmetric five-peaked toneburst signal generated is defined where the minimum value 
occurs after the maximum (Figure 7). and were checked to ensure that they occur on 

subsequent peaks by Eq. (2). If and did not satisfy Eq. (2), a secondary algorithm 

must be implemented to determine • 

Lin ~ Lax < J ( 2 ) 


where N equals the number of peaks (5) and is the center frequency (300 kHz). 




For each of the twelve measurements, and K(s)\^ did satisfy Eq. (2). was then 

calculated as the average of and (Eq. (3)). 

^a(s) - ^ tj; 


A plot of the actuation waveform and one of the sensor responses with the key time points in the 
algorithm is shown in Figure 7 and Figure 8, respectively. 
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Figure 7. Actuation waveform. Figure 8. Plate response from sensor located 

five inches from actuator. 

The average group velocity was then determined from the group velocities of each of the twelve 
actuator-sensor pairs. The three points in Figure 9 correspond to the average group velocity of 
each of the three panels. The cure time value is denoted by the time at which the composite panel 
began to cool during the curing cycle (Figure 1). The error bars show the variation based on a 
90% confidence interval of the twelve measurements. The average group velocity increased by 
5.52 percent from Panel #1 to Panel #2 and another 1.26 percent from Panel #2 to Panel #3 (Figure 

9). 



Cure Time (mins) 

Figure 9. Group velocity normal to the fiber direction on each of three panels of varying final 

degrees of cure. 

It is expected that each subsequent panel with longer cure times leads to a higher final degree of 
cure and lower porosity level. As shown in Figure 10, a decrease in the number of voids was 


observed from each subsequent panel especially between Panel #1 and Panel #2. A lower porosity 
level in each subsequent panel can be confirmed from micrographs and/or acid digestion which 
were not completed at this time so that the panels could be utilized for future tests. 



Panel #1 Panel #2 Panel #3 

Figure 10. C-scan images of fabricated panels. 

The average group velocity increased monotonically from Panel #1 to Panel #3. Therefore, 
preliminary results indicate that the group velocity of a guided wave could possibly be used as an 
initial index (or metric) to determine the degree of cure, and infer the level of porosity. This 
process can be transitioned into an in-line detection system that will be able to monitor the degree 
of cure and porosity in real-time within a part. The system will consist of high temperature 
actuators and sensors that can withstand the temperature and pressure during cure and high 
temperature cables that will connect the actuators and sensors to the electronics outside the cure 
environment. Many technical details will need to be determined to make this transition including 
mounting of the actuators and sensors to the mold, coupling the energy from the sensor to the part, 
and filtering out the effect of temperature on the signals. 

3.2 Overlap Detection 

After the panel was cured it was removed from the mold and C-scanned with the same system 
described in Section 2. 1 to establish a baseline to compare with the guided wave approach (Figure 
11 ). 


Overlap 



Figure 11. C-scan image of panel with intentionally introduced overlap and gap. 

For the guided wave approach, the area that is to be imaged (46 mm x 26 mm) to reveal the overlap 
region in Figure 6 depends on the aperture of the air-coupled transducer. The angle of incidence 
of the air-coupled transducer was set to 13.5 degrees from normal. This value is optimized to 
ensure maximum amplitude of the desired mode of the propagating wave in the plate [22] and can 
be determined from the phase velocity dispersion curve (Figure 12). 



Figure 12. Phase velocity dispersion curve for the [O 4 / 9 O/ 9 O 2 / 9 O/O 4 ] composite plate. 

At the center frequency of the excited wave (250 kHz), two fundamental wave modes 
simultaneously exist (Ao and So). Ao is the antisymmetric and So is the symmetric zero-order Lamb 
wave mode. Since the magnitude of the transverse velocity generated by the mode is much 
smaller than the Ao mode, the 5'o mode captured by the LDV is much less sensitive. Therefore in 
this study only the Ao mode was chosen for processing the wave signals. At 250 kHz, the phase 
velocity of the Ao mode is 1.452 km/s. Using Eq. (4) [22], the optimal angle of excitation, 0 , can 
be calculated 


(2 

sin0 = — 


( 4 ) 


where c is the sound velocity in the excitation medium (Air: 343 m/s) and is the phase velocity 

of the Aq mode in the plate (1.452 km/s). The resulting optimal angle of excitation for the air- 
coupled transducer is 13.66 degrees (13.5 degrees was used in the experiment). 

The out-of-plane velocity of the wavefield, defined by w{x,y,t), can be assembled from all the time 
series from the LDV in the scan area. The wavefield is then windowed in time after the actuation 
is triggered, and prior to the wave reflected back from right edge of the plate. The wavefield 
w(x,y,t) is plotted at several time steps as a 3-D surface plot (Figure 13) and as a 2-D contour plot 
(Figure 14). 



Figure 13. 3-D surface plots of the out-of-plane velocity, w(x,y,t), of wavefield at times (a) 155.9 
p5 , (b) 165.6 p 5 , (c) 175.4 pj' , (d) 185.2 pj' , (e) 194.9 ps , and (f) 204.9 pj' . 














Figure 14. 2-D contour plots of the out-of-plane velocity, w{x,y,t), of wavefield at times (a) 155.9 
, (b) 165.6 P5' , (c) 175.4 , (d) 185.2 p.? , (e) 194.9 p^' , and (f) 204.9 p.s in the scan area 

where the red rectangle outlines the imaged area shown in Figure 6. 

The cumulative total wave energy (CTWE) can be calculated by Eq. (5) [19] 

CTWE{x,y) = W^{x,y,t) At 


( 5 ) 


where t. is the time at the beginning of the time window measurement and is the time at the 
end of the time window (prior to the back reflection from right edge of the plate). 

CTWE signifies the area illuminated by the air-coupled transducer and serves as a comparison to 
the area imaged by the ZLCC imaging condition. 


Figure 15. Cumulative total wave energy (CTWE) in the imaged area (outlined by red rectangle 

shown in Figure 6). 

The ZLCC imaging condition for determining the size of the defect is based on the concept that 
the incident and backscattered wave are in phase on the damage boundaries [20]. From the 
experimental setup (Figure 5), it can be observed that the incident wave propagates to the right 
(forward propagating) and the backscattered wave propagates to the left (backward propagating). 
To decompose w{x,y,t) into the incident and backscatter wave, w{x,y,t) is first transformed from 
the time-space domain to the frequency- wavenumber domain [19]. This can be accomplished 
using a two-dimensional Fourier transform (Eq. (6)). 


Note the sign convention where a positive denotes a wave propagating forward to the positive 
(or right) x-direction whereas a negative denotes a wave propagating (scattering back) to the 
negative (or left) x-direction. 

The total wavefield can be readily separated in the frequency-wavenumber domain into the 
incident W.{k^,y,(a) and the backscatter W^^(k^,y,a>) wave. Both the incident wave and the 

backscattered wave can be inverted to the frequency-space domain using a one-dimensional spatial 
Fourier transform (Eq. (7)). 
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For each measurement location (x.,y.), the ZLCC value, is found using Eq. (8) and Eq. (9) 

[ 20 ] 



( 8 ) 

( 9 ) 


I{x.,yj) = Yjy(.x.,y.,(Hj^) 

CO 


where w^^(x.,y.,ca) is the complex conjugate of W^^(x.,yj,co). 

Figure 16 shows the ZLCC value for the imaged area with measurements in a 1 mm x 1 mm grid. 
The overlap location is denoted by a filled red rectangle below. 
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Figure 16. Zero-lag cross-correlation in the imaged area from Figure 6 (overlap location is 

denoted by the filled red rectangle). 

For both Figure 17 and Figure 18, /(x, , is calculated with white noise removal at each step 
of /(x;,v^.,(o^) before summation. The white noise was mitigated by using a two-dimensional 

single wavelet transform (Daubechies 10 level 5). The horizontal, vertical, and diagonal details 
coefficients were calculated using Matlab’s Wavelet Toolbox. Figure 17 shows the result when 
the measurements were taken in a 1 mm x 1 rnm grid. Figure 18 shows the result when the 
measurements were taken in a 2 mm x 2 mm grid. In both figures, the overlap location is denoted 
by a filled red rectangle. 




Figure 17. Zero-lag cross-correlation domain with white noise removal at each frequency step 
before summation in the imaged area from Figure 6 (overlap location is denoted by the filled red 

rectangle). 
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Figure 18. Zero-lag cross-correlation domain with white noise removal at each frequency step 
before summation in the imaged area from Figure 6 (overlap location is denoted by the filled red 

rectangle). 

The ZLCC imaging condition was utilized to estimate the overlap region within the imaged area 
up to a 2 mm x 2 mm grid size. Using a two-dimensional single wavelet transform (Daubechies 
10 level 5) to mitigate the noise from the image leads to a sharpened image of the overlap location 
and reduces the presence of false-positives. From Figure 17 and Figure 18, one can notice that the 
imaged region was biased toward the right side of the overlap region. A possible reason for this 
outcome is that the wave is traveling from left to right with all of its energy in the original twelve 
layers. At the left boundary, two prepreg layers are added in the ninety-degree direction 
(perpendicular to wave propagation direction) causing some degree of the waves reflected (or 
scattered). However, higher backscattered waves are induced at the right boundary of the overlap 
where two of the original prepreg layers are suddenly dropped. 




The accuracy of the result is limited by following factors. First, the actuation frequency is limited 
by the response sensitivity of the LDV. If a higher frequency could be measured, then the 
propagating wave would have a smaller wavelength leading to a greater interaction with the 
overlap. Second, only a single excitation source location was used. The excitation source could 
be moved up/down the plate surface and a similar imaged area scanned. The results could then be 
superimposed to further highlight the overlap region. Finally, the Ao mode is a flexural mode; 
therefore, the stress is at a minimum at the neutral axis of the plate. This causes the wave 
interaction with the overlap to be minimized when the overlap location is at the center of the plate 
compared to any other location through the thickness. 

4. SUMMARY OF RESULTS AND CONCLUSIONS 

Three 24-ply unidirectional composite panels were fabricated using a vacuum press where each 
subsequent panel had a higher final degree of cure and expected lower porosity level than the 
previous one by varying the cure cycle. Guided waves were propagated and recorded in a pitch- 
catch configuration normal to the fiber direction using piezoelectric actuators and sensors (Figure 
2 and Figure 3). From the group velocity calculated from the actuator-sensor pairs, the average 
group velocity was determined. The average group velocity increased from 1 .53 km/s to 1 .62 km/s 
to 1 .64 km/s from Panel #1 to Panel #2 to Panel #3 (Figure 9). Thus, the group velocity of a guided 
wave may have the potential to serve as an initial index for estimating the degree of cure and level 
of porosity with further work including micrographs and acid digestion. 

Ultimately, the future work will be to fully transition to an in-line detection system that will be 
able to monitor the degree of cure and porosity in real-time within a part based on group velocity. 
The in-line detection system will utilize high temperature actuators and sensors for guided wave 
emission and response detection as well as investigate other metrics such as attenuation of wave 
energy which will play a critical role during the early stages of cure. 

A composite panel with an intentionally introduced overlap in the middle two layers was cured in 
a vacuum press. A stationary unfocused air-coupled transducer excited a guided wave in the 
composite plate and a LDV recorded the velocity of the plate in the transverse direction. The 
response was transformed from the time-space domain to the frequency-wavenumber domain and 
then separated into the incident and backscatter wave. A zero-lag cross-correlation (ZLCC) 
imaging condition was proposed to image the overlap region where the incident and backscattered 
wave are in phase (the boundaries of the overlap). A two-dimensional single wavelet transform 
was used to mitigate the noise from the ZLCC matrix at each frequency step before summation. 
The location of the overlap was identified within the imaged area with a grid size up to a 2 mm x 
2 mm (Figure 17 and Figure 18). 

The future work concerning this detection system will be to mitigate/modify the factors identified 
in Section 3.2 to sharpen the image within the constraints of the equipment and extend the 
technique for the detection and size determination of other composite manufacturing defects 
including gaps, discrete locations of high porosity, voids, through-the-thickness fiber waviness, 
and in-plane fiber waviness. 
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